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Pion decay constant in quenched QCD with Kogut-Susskind quarks* 
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We report results on the flavor breaking effect in the pion decay constant f„ calculated with the Kogut-Susskind 
quark action. Numerical simulations are carried out at j3 = 6.0 and 6.2. We find that the use of non-perturbative 
renormalization factor leads to results in the continuum limit well convergent among various KS flavors. As the 
best value for f-^ we obtain = 89(6) MeV from the conserved axial vector current. 



1. Introduction 

The Kogut-Susskind (KS) quark action has the 
well-known feature that SU(4) flavor symmetry 
is broken down to U(l) subgroup at finite lattice 
spacing. The restoration of full flavor symmetry 
toward the continuum limit has been previously 
examined for pion mass jl],|2|] . Here we extend the 
examination to the pion decay constant by com- 
paring results for various KS flavors to that in 
the U(l) channel for which the renormalization 
constant equals unity. 

This 

comparison is made both for perturbative H and 
non-perturbative renormalization factors, the lat- 
ter evaluated Q with the method of Ref. ||. 

Numerical simulations are carried out in 
quenched QCD at (3 = 6.0 and 6.2 employing 
32 3 x 64 and 48 3 x 64 lattices. Other lattice pa- 
rameters are summarized in Tabic []]. 

2. Formalism and simulation 

In the hypercubic notation for the KS fermion 
fields, the axial vector current in the KS flavor 



Table 1 

Lattice parameters of simulation. 
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L A 


T 


m q a 


#conf. 


a- x (GeV) 




6.0 


32 3 


64 


.010-.030 


100 


1.93(2) 


.53-.72 


6.2 


48 3 


64 


.008-.023 


60 


2.70(5) 


.55-. 75 



channel F has the form 
K = ® &)<!>■ 



(1) 



The fields 4> and (f> in ([l]) are generally separated 
over a 4-dimensional hypercube. We use both 
gauge-invariant and non-invariant operators, in- 
serting the product of link variables between <j> 
and 4> for the former, and employing the Landau 
gauge fixing to deal with the latter. 

The one-loop perturbative results for the renor- 
malization constants have been worked out 
m Ref. @. We use tadpole-improved values em- 
ploying the tadpole-improved MS coupling with 
q* = 1/a; they will be denoted as Z A 

(N)F 



For non 

perturbative renormalization constants Z A 
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we 

take results of Ref. ]g at {apf = 1.0024." 

Quark propagators are calculated for 16 wall 
sources, each corresponding to a corner of a hy- 
percube. These propagators are combined to 
form the correlator (A^.{t)-Kyy{fy) f° r each flavor 
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Figure 1. Continuum limit of pion masses. 

F, where ir^y is the pion field for the wall source, 
thereby enhancing signals 

We find the quark mass dependence of (m^) 2 , 
f£ and vector meson masses to be well described 
by a linear function. Hence a linear fit is employed 
for the chiral extrapolation in the quark mass m q . 
The scale is set by the p meson mass in the spin- 
flavor channel <8> £fc- Errors are estimated by 
the single elimination jackknife procedure. 

3. Pion masses 

We show in Fig. |l| values of (m^) 2 at m q — 
as a function of a 2 . The 16 KS flavors 
are grouped into 8 irreducible representations 
(irreps) given by@ (1-d irreps) 

and 6fe£5,6fe£4,&£m,£fc (3-d irreps). We ob- 
serve that the 8 irreps form a degeneracy pat- 
tern f s , (fkfsjkfs), &£m), (£4, €k),I at fi- 
nite lattice spacing. This feature was initially 
observed numerically in Ref . H , and recently the- 
oretically explained in Ref. M. We also see that 



non-zero values of 



,F\2 



for the channels other 



than £ 5 vanish quadratically in a toward the con- 
tinuum limit as expected. 

4. Pion decay constants 

In Fig. H we illustrate how the bare values of 
f£ change under renormalization, taking results 
for gauge invariant current at (3 = 6.2. The bare 
values for the 8 irreps again form a degeneracy 
pattern. The pattern reflects the distance of <fi 
and (f> in ([!]), and is different from that of m^. 
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Figure 2. obtained with gauge-invariant ax- 
ial vector current as functions of m q at (3 = 6.2 
before renormalization (top), with perturbative 
(middle) and non-perturbative (bottom) renor- 
malization factors. 

The perturbative renormalization constants 
(middle frame in Fig. |^) help to reduce the dis- 
crepancy among the bare values of f£ . It is clear, 
however, that the discrepancy is much more re- 
duced with the non-perturbative renormalization 
constants (bottom frame). 

In Fig. [3| we plot the continuum extrapola- 
tion of the pion decay constants obtained with 
gauge-invariant currents. The top figure shows 
results with perturbative renormalization factors, 
and the bottom figure with non-perturbative fac- 
tors. Lines are quadratic fits in a 2 following the 
0(a 2 ) scaling violation expected for the KS quark 
action. Similar figures for the decay constants ob- 
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Figure 3. Continuum limit of f£ with perturba- 
tive (top) and non-perturbative (bottom) renor- 
malization factor obtained with gauge-invariant 
axial vector current. 

tained with the gauge non-invariant currents are 
shown in Fig. [|. 

The gauge-invariant current for the flavor £5 
does not require renormalization, from which we 
obtain / w = 89(6) MeV in the continuum limit. 
This result is consistent with the experimental 
value of 92.4(3) MeV||; possible deviations due 
to quenching is not visible within our 7% error. 

We observe in Figs. || and [| that the pertur- 
bative renormalization factors at one-loop order 
are not sufficient to ensure restoration of fla- 
vor symmetry in the continuum limit. In con- 
trast, the decay constants evaluated with the non- 
perturbative renormalization factors agree much 
better already at finite lattice spacings, and their 
continuum limits are well convergent, the central 
values coinciding with each other well within the 
errors. 
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Figure 4. Same as Fig. g obtained with gauge 
11011- invariant axial vector current. 
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